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ABSTRACT 

The quadrupole power of cosmic microwave background (CMB) temperature anisotropies seen in the 
WMAP data is puzzHngly low. In this paper we demonstrate that Minimum Variance Optimization 
(MVO), a technique used by many authors (including the WMAP science team) to separate the 
CMB from contaminating foregrounds, has the effect of forcing the extracted CMB map to have zero 
statistical correlation with the foreground emission. Over an ensemble of universes the true CMB 
and foreground are indeed expected to be uncorrelated, but any particular sky pattern (such as the 
one we happen to observe) will generate non-zero measured correlations simply by chance. We call 
this effect "cosmic covariance" and it is a possible source of bias in the CMB maps cleaned using the 
MVO technique. We show that the presence of cosmic covariance is expected to artificially suppress 
the variance of the Internal Linear Combination (ILC) map obtained via MVO. It also propagates 
into the multipole expansion of the ILC map, generating a quadrupole deficit with more than 90% 
confidence. Since we do not know the CMB and the foregrounds a priori, there is therefore an 
unknown contribution to the uncertainty in the measured quadrupole power, over and above the 
usual cosmic variance contribution. Using the MVO on a series of Monte Carlo simulations that 
assume Gaussian CMB fluctuations, we estimate that the real quadrupole power of the CMB lies in 
the range [305.16, 400.40] (^fK^) (at the 1 - ct level). 

Subject headings: cosmology: cosmic microwave background — cosmology: observations — methods: 
data analysis 



1. INTRODUCTION 

Detailed measurements of the pattern of temperature 
anisotropies in the cosmic microwave background (CMB) 
have provided cosmologists with unprecedented oppor- 
tunities to probe the physics of the early Universe. In 
particular, the Wilkinson Microwave Anisotropy Probe 
(WMAP) data (Bennett et al. 2003b,c; Hinshaw et al. 
2003b, 2007) have played a pivotal role in the establish- 
ment of the standard "concordance" cosmological model 
and have heavily constrained alternatives to the infla- 
tionary paradigm for the origin of large-scale structure. 

According to the concordance cosmology, the varia- 
tions in temperature of the CMB across the celestial 
sphere should possess a very broad spectrum usually ex- 
pressed in terms of the amplitude of spherical harmonic 
modes labelled by the usual angular harmonic frequency 
£. The power of the pattern of anisotropies power is ex- 
pected to be strongest for the quadrupole {£ — 2) which 
correspond to variations on an angular scale of 90° on the 
sky. In principle, therefore, the quadrupole should be the 
easiest harmonic mode to detect. In practice, however, 
the measured quadrupole anisotropy of the CMB from 
WMAP seems to sit rather uncomfortably with the stan- 
dard cosmology: its amplitude seems too low (Efstathiou 
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2003a; de Oliveira-Costa et al. 2004). Of course, the 
power that resides in each multipole of the CMB pattern 
(including the quadrupole) is not so much measured as 
extracted. There are various diffuse foreground emissions 
to be separated from the "true" primordial CMB, such 
as free-free emission arising from electron-ion scattering, 
synchrotron emission from cosmic ray electrons acceler- 
ating in the galactic magnetic fields, and dust emission. 
However, although a number of authors have employed 
different foreground subtraction schemes, all find the de- 
rived quadrupole power to be lower than the Concor- 
dance Cosmological Model with a significance level of 
5% or better (Efstathiou 2003a; Eriksen et al. 2004; de 
Oliveira-Costa & Tegmark 2006; Park, Park & Gott 2007; 
Saha et al. 2007). This has even prompted suggestions 
of physics beyond the standard model (Efstathiou 2003b; 
Luminet et al. 2003). 

In this paper we suggest that the quadruple power 
deficit in WMAP may be an artifact of the Minimum 
Variance Optimization method used by many authors to 
clean CMB maps from foreground contamination. In the 
next Section, we explain the basics of the MVO method. 
In Section 3 we show how this method is prone to a bias 
introduced by the accidental alignment of structures in 
the CMB and foreground templates leading to an effect 
we call "cosmic covariance" . In Section 4 we quantify 
the likely effect of this bias using Monte Carlo simula- 
tions, and we estimate the quadrupole power taking into 
account the effect of cosmic covariance. We present our 
conclusions in Section 6. 

2. MINIMUM VARIANCE OPTIMIZATION 

Although various schemes have been proposed for sub- 
tracting foregrounds from CMB data, one concept that 
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tends to be in common for multi-frequency cleaning 
methods is Minimum Variance Optimization (MVO). 
The theoretical basis for MVO is the primordial CMB 
temperature anisotropies arise from black-body radiation 
and therefore constitute a frequency-independent signal 
that persists across different observed frequency bands 
(Tcgmark & Efstathiou 1996). The observed sky tem- 
perature in a given band i is therefore just 

= To + (1) 

where is the CMB signal and fi the foreground emis- 
sion band i, which will bo a composite of dust, free free 
and synchrotron sources. For the purposes of this pa- 
per we assume that instrument noise is negligible. The 
frequency-independent signal can be flushed out by an 
Internal Linear Combination (ILC) from the maps of M 
frequency bands with weighting coefficients Wi, where 
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1 (Bennett et al. 2003c): 
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iTiip) = n{p)+y^wji{p) = Te(p)+Tfr(p), 

(2) 

where T^c {p) denotes the ILC map value at a pixel p and 
Tfr the foreground residual. In general, 

Var(Tue) = Var(rc) + Var(Tfr) + 2Cov(re, Tf,), (3) 

where the variances and covariance are measured over 
an ensemble of possible skies. The MVO technique is 
based on the a priori assumption that the frequency- 
independent signal (i.e. the true CMB) should be sta- 
tistically independent of any foreground contribution, so 
that the covariance term vanishes and 



Var(Tiic) = Var(Tc) + Var(Tfr), 
or in other words, 
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where the angle brackets denote ensemble averages as 
above. Minimizing the variance of the ILC map is then 
equivalent to minimizing the foreground contamination 
of a map of thermal noise. 

The WMAP 3-year ILC map and, more importantly, 
the power spectrum for £ < 30 are"'' thus produced by 
employing the MVO on 5 frequency maps K (23 GHz), 
Ka (33 GHz), Q (41 GHz), V (61 GHz) and W (94 GHz) 
bands in the pixel domain for 12 separate regions, where 
Region marks the largest region with \b\ > 15° (about 
89% of the whole sky), and Regions 1 — 11 are those 
around the Galactic plane. For each region the 5 fre- 
quency band maps are linearly combined and a set of 

weighting coefficients for each region are obtained 
in such a way that the resultant variance is forced to be 
minimum: 
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Note that Equation (5) involves the requirement of sta- 
tistical independence over an ensemble of universes. As 

^ For power spectrum for £ > 30, WMAP science team applies 
a foreground template fitting method. 



we have only one universe (i.e. one CMB sky), the calcu- 
lation of variances and covariances can only be done over 
the set of pixels corresponding to a single sky, and not 
over an entire probably distribution. The sky covariance 
between Tc and Tfr need not be exactly zero, since there 
is always a chance of some chance coincidence of features 
such as hotspots between foreground and true CMB. If 
this happens then 
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As we shall now show, this leads to a bias in the recovery 
of the CMB map. 

3. COSMIC COVARIANCE 

Let us assume (for illustrative purposes only) that the 
components of the foreground emission at different fre- 
quencies have the same spatial distribution (i.e uniform 
foreground spectra). In reality the foregrounds at dif- 
ferent frequency bands have spatial variations of the 
spectral indices, which can amplify the resulting bias. 
In the uniform case we have fi{p) = SiF{p), where 
Si = S{vi) is the composite frequency spectrum of the 
foreground emission and F is the common spatial distri- 
bution over the set of pixels. Thus frequency band maps 
Ti{p) = Tc{p) + SiF{p) and the ILC map 



Tnc{p)=T,{p)+TF{p), 



(8) 



where F = WiSi. The variance of the ILC map can 
be written as follows (Hinshaw et al. 2007) 



ilc - 

where 
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a,,F = {TcF)-{n){F) 



(10) 

is the "cosmic covariance" between the CMB and the 
foreground spatial distribution F. In this and the subse- 
quent equations, the angle brackets denote averages over 
the pixels of a single sky rather than ensemble averages; 
c.f. Equation (5). 

The term "cosmic covariance" is coined similarly to 
"cosmic variance" (White, Krauss and Silk 1993) because 
both effects emanate from the fact that we must do our 
statistical analysis using a single version of the celestial 
sphere. Cosmic variance arises from the fact that a sin- 
gle sky does not offer sufficient spherical harmonic modes 
to make an accurate estimate of the ensemble-averaged 
power at low multipoles. In the case of cosmic covariance, 
the single sky prevents us from being sure that our calcu- 
lation of the correlation between foreground and CMB is 
exact, particularly for the quadratic minimization. The 
value (7c, F will be non-zero merely because the CMB and 
the foregrounds happen to line up to a certain degree 
simply by chance. If both CMB and foreground contain 
large-scale structures then the effect of cosmic covariance 
will be larger as the sky will contain fewer independent 
regions. 

The cosmic covariance cTc^f can also be expressed via 
the correlation coefficient Xcf- ctcF = XcfCqO'f, where 
XcF, the "cosmic correlation coefficient" characterizes 
the level of resemblance in morphology between the CMB 
and the foreground spatial distribution, and — 1 < X^f < 
1. 
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TABLE 1 

Correlation coefficients X between the full-sky WMAP 
forec;round maps (lower left triangle) and between those 

IN THE WMAP -defined REGION 0, WHICH ACCOUNTS FOR 89% 
OF THE FULL SKY (UPPER RIGHT TRIANGLE). OnE CAN SEE THAT 
K, Ka, Q AND V CHANNELS (EXCEPT W) HAVE CORRELATION 
VERY CLOSE TO UNITY. AlL OTHER REGIONS OF THE WMAP 
CHANNELS EXCEPT W SHOW THE SAME TREND. 

Employing the MVO criterion to the variance of the 
ILC map, we get 
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-XcfCc/cf and the ILC map is then 
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One can see if cosmic covariance is zero (so is Xcf), then 
Tnc{p) = T,{p). 

The last term —XcfF{p)(7c/'^f is usually discussed in 
terms of an anti-correlation bias (Hinshaw et al. 2007), 
but it is in fact due to the MVO serving to eliminate 
the covariance between the ILC and the foreground spa- 
tial distribution. To see this, note that the covariance 
between the ILC and the foreground is forced to vanish: 

^iicF = (TilcF) - {Tnc){F) = ac,F - = 0. (13) 

ctf 

The ILC map is thus guaranteed to have a morphology 
that has no resemblance to that of the foregrounds, re- 
gardless of whether or not the CMB has any such resem- 
blance. That is to say, even if there exists any cosmic 
covariance between the true CMB signal and the fore- 
grounds, it is duly subtracted by the MVO, and the re- 
sulting ILC map must display zero sky correlation with 
the foregrounds. 

Forcible subtraction of the unknown cosmic covariance 
also affects the overall variance. Putting F = — Xcfctc/ctf 
into Eq.(9), we have 

al=a',{l-X,F'). (14) 

Subtraction of the cosmic covariance manifests itself as a 
deficit in the total variance of a factor (1 — Xcf^). This 
can partially account for the low variance in the CMB 
map (Monteserin et al. 2007). As we do not know Tc 
and F a priori, we are left with no information about 
the level of cosmic covariance after it is subtracted. 

One of the dangers of the MVO approach is that it 
encourages circular reasoning. The a priori assumption 
is that the CMB and the foregrounds are statistically 
independent. The MVO produces an ILC map that has 
zero sky correlation with the foregrounds. This seems to 
confirm the correctness of this assumption. However, the 
MVO must produce this result whatever the morphology 
of the input templates. 

Although the above analysis is based on the existence 
of composite foregrounds with uniform frequency spec- 
tra, we can see that it is not far from reality. In Ta- 
ble 1, we show the pixel- by- pixel correlation coefficients 
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Fig. 1. — Decreasing of correlation with the foreground due to 
the MVO. We show the histogram of X^p and X^p, the correlation 
coefficient between the foreground and the simulated CMB (left) 
and the retrieved CMB (right) via the MVO. 

between the WMAP composite foregroimds in WMAP- 
defined Region 0, which contains ~ 89% of the whole 
sky (upper right triangle). One can see that all channels 
have excellent agreement in terms of morphology (with 
the exception of W channel due to synchrotron emission) . 
The lower left triangle part shows correlations between 
full-sky foreground maps, which also display the same 
trend. 

4. MONTE CARLO SIMULATIONS 

To demonstrate the power deficit and reduction of cos- 
mic covariance in the ILC map due to the MVO, we 
perform simulations of 2000 sets of frequency band maps 
corresponding to the WMAP K, Ka, Q, V and W chan- 
nels. These comprise the sum of a simulated CMB sig- 
nal (assuming Gaussian fluctuations), mock foreground 
maps and instrument noise maps. The CMB signal is 
simulated with the WMAP best-fit ACDM model. The 
mock foreground maps at WMAP 5 frequency bands are 
produced as follows: we take the WMAP composite fore- 
ground at K band, which is multiplied by a factor pi such 
that piffK = o'iji = Ka, Q, V and W where ai is the stan- 
dard deviation of the WMAP composite foreground maps 
(in thermodynamic temperature units). As such, these 
mock foreground maps have exactly the same morphol- 
ogy, but each has the same variance as the WMAP fore- 
ground map of the corresponding frequency band. The 
Gaussian instrument noise map is with = 2 x 10~^ 
(mK^). To follow the WMAP procedure, both the simu- 
lated CMB and noise maps are smoothed to 1° FWHM 
before summation with the mock foregrounds (which are 
already smoothed). The MVO is then employed on the 
internal linear combination of the band maps to retrieve 
the CMB signal. 

In Figure 1 we plot the histogram of Xcf (left), the 
correlation coefficient between the simulated CMB and 
the common foreground map, and that of VjF (right), the 
correlation coefficient between the retrieved CMB (ILC) 
and the common foreground map. One can see the sky 
corr(^lation is reduced below a level 2 x 10~* after the 
appfication of the MVO. 

The deficit in the variance XcF^a"^ is registered in the 
angular power spectrum Q via af^^ = (47r)~^ S^(2^ + 
l)Cf'^. Note that Ce is always positive, so the deficit in 
variance manifests itself as a few dips compared to the de- 
sired CMB angular power spectrum Cg. Since the MVO 
acts on the overall variance, it only ensures the variance 
is minimum; the result for individual spherical harmonic 
modes is less obvious. Decomposing Eq.(12) into spher- 
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Fig. 2. — Normalized histograms (thin) and cumulative his- 
tograms (thick) of the power difference A(<5T|) = STf^ — STf ^ 
from 2000 realizations for ^ = 2, 3, 4 and 5, where ST"^^ denotes 
the simulated CMB power at i mode and STf. the retrieved CMB 
(ILC) power. The MVO is employed on frequency band maps with 
the foregrounds having the same morphology (see text for the sim- 
ulations). One can see that the CMB retrieved via MVO has a 
skew distribution for t = 2, indicating the power of the retrieved 
CMB is smaller than the simulated ones, 
can examine its angular power spectrum 
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One can see that it is possible for some multipoles to 

acquire excess power (bumps) over if, for instance, 
tt/2 < (pj^ — < 2,11/2. So in order to preserve the 
overall deficit, it shall drag the dips at other multipoles 
even further down. 

In Figure 2 we plot the normalized histogram and cu- 
mulative histogram of the power difference l^{5Tf) = 
6Tl^-STl^, where ST^ = £(€+l)Q/(27r) is the tempera- 
ture anisotropy power at multipole £, the subscripts s and 
i denote the simulated CMB and the retrieved CMB from 
MVO, respectively. It is easy to see that, particularly for 
multipole number 1 = 2, the simulation ensemble forms 
a skewed distribution with a tail on the negative side, in- 
dicating the power of the retrieved CMB is more likely to 
be smaller than that of the input (simulated) CMB. For 
other multipoles the skew is less obvious. This shows, 
therefore, that the deficit in variance due to the MVO is 
most likely registered in the quadrupole. 

One should note that the power deficit claimed by the 
WMAP team is calculated from Region (Full sky with 
Galaxy cut) and then corrected by using a maximum 
likelihood estimate (Hinshaw et al. 2007). The power we 
have calculated above, however, is directly from full-sky 
simulations because one can then single out the effect 
of the MVO, whereas extrapolation of the power from 
incomplete sky coverage may introduce yet another sys- 
tematic error, particularly for large angular scales. 

5. ESTIMATION OF THE QUADRUPOLE POWER 

Since we know the cosmic covariance is going to be 
subtracted by the MVO, there is an unknown error in 
the quadrupole power even before the cosmic variance 
becomes involved in the interpretation of the data. Con- 
sequently, the quadrupole power we can estimate is at 
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Fig. 3. — Distribution of the CMB quadrupole power. We employ 
the MVO on frequency band maps of the WMAP K, Ka, Q, V 
channels (based on the excellent agreement in morphology between 
the foregrounds) and add on the estimated estimated quadrupole 
the missing power from 10'' realizations of 4-cliannel MVO. The 
true CMB quadrupole falls in [305.16, 400.40](mK2) (at the 1 - <t 
level) . 

best not going to be a single value, but a sampling dis- 
tribution which we construct with help of Monte Carlo 

simulations, assuming the CMB is a Gaussian random 
field. According to Table 1 (lower left triangle), the full- 
sky WMAP foregrounds at K, Ka, Q and V channels 
(except W) have correlation coefficients deviated from 
unity by less than 5 x 10~^. We can employ the MVO 
on these 4 channels to retrieve the quadrupole power and 
then add on the missing power from Monte Carlo simu- 
lations. 

The weighting coefficients for the linear combina- 
tion from the 4 frequency band maps are (6.2419 x 
10~-\ -0.21597, -0.37970, 1.5894) and we denote the re- 
sultant map as ILC4. The correlation coefficients be- 
tween the ILC4 and the derived foregrounds (frequency 
band maps subtracting ILC4) at K, Ka, Q and V are 
X = (8.3553,8.3510,8.3351,8.2826) x 10"^, respectively, 
which confirms the presence of negligible correlations af- 
ter the application of MVO. The quadrupole power from 
the ILC4 map is ST^ = 326.30 /iK^. We then add the 
missing power for quadrupole from 10'' realizations of ex- 
actly the same procedure as described earlier in this Sec- 
tion, except this time using WMAP Q band foreground 
as common foreground for the 4 channels. In Fig. 3 we 
plot the histogram of the resulting estimates of true CMB 
quadrupole power, which falls in [305.16, 400. 40](/iK2) 
(at the 1 — a level) . 

6. DISCUSSION AND CONCLUSIONS 

We have shown that the MVO method for eliminat- 
ing foreground contamination from CMB maps is prone 
to a bias if there exist accidental alignments between 
foreground features and structures in the true CMB sky. 
This bias acts in such a way that a suppression of the 
quadrupole is more likely than an enhancement. It also 
reduces the overall variance of the recovered map com- 
pared to the true CMB sky. This suggests that attempts 
to attribute the low quadrupole to exotic physics may be 
premature. 

The MVO was originally designed for the separation 
of Gaussian signals, whose statistical properties are com- 
pletely characterized by their power spectra. For a linear 
combination of such signals, the variance is clearly the 
natural choice for a parameter to be minimized, and in 
such a case this method generally leads to a good recon- 
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struction of the power in each component. In this pa- 
per, we have demonstrated even for the MVO to work 
on Gaussian signals it requires a statistical ensemble: 
for a single realization the role of cosmic covariance be- 
comes crucial. The foregrounds relevant to CMB analysis 
are also highly non-Gaussian, which makes this problem 
even worse. 

For composite foregrounds with varying spectra, the 
variance after application of the MVO will induce even 
more foreground contamination in the ILC map because 
the weighting coefScients will then be tuned not only 

for the subtraction of cosmic covariance, but for vary- 
ing morphology among the different foregrounds. It is 

® http://healpix.jpl.nasa.gov/ 
http : //wn . glesp . nbi . dk/ 



then possible that the variance deficit generated by the 
subtraction of cosmic covariance can be compensated by 
a foreground residual, which can push the quadrupole 
power back to the value that is consistent with the Con- 
cordance Cosmological model. Extra caution, therefore, 
will have to be applied when employing the MVO to the 
interpretation of results from the upcoming ESA Planck 
Experiment. 
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